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We construct a qubit algebra from field creation and annihilation operators acting on a global 
vacuum state. Particles to be used as qubits are created from the vacuum by a near-deterministic 
single particle source. Our formulation makes the space-time dependence of the qubits explicit, 
preparing the way for quantum computation within a field framework. The method can be general- 
ized to deal with interacting qubits whose wavepackets are not perfectly matched to each other. We 
give an example of how to calculate the Heisenberg evolution of a simple two-qubit circuit, taking 
expectation values in the field vacuum state. 
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I. INTRODUCTION 

Much of quantum information theory is couched in the 
formalism and language of the Schrodinger Picture of 
quantum mechanics, in which the qubit states explicitly 
evolve as they pass through circuits or down communi- 
cation channels pQ . Key properties such as entanglement 
depend specifically on the Hilbert space structure of the 
qubit state after processing. In contrast, a completely 
equivalent approach is to work in the quantum mechani- 
cal Heisenberg Picture, in which it is the Pauli operators 
describing observables that evolve, whilst the state of the 
system is always the global ground-state of the qubits [5] . 
In this formalism the properties of the system depend on 
the Pauli operator algebra and specifically on its commu- 
tator structure. Heisenberg evolution is used extensively 
in the stabilizer formalism for error correction and graph 
states [3]. 

A limitation of current quantum information theory is 
that it treats qubits as stationary point particles possess- 
ing no space-time extent or dynamics. A specific excep- 
tion is in optics where the finite extent and propagation 
of photon modes has been included in linear optical quan- 
tum computation circuits [1]. 

In this paper we present a recipe for performing quan- 
tum information calculations in the Heisenberg Picture 
in which the Pauli operators are explicitly constructed 
from field operators with a spacetime structure and the 
initial state is always the global field ground, or vacuum 
state. This is in contrast to the usual approach, where 
the ground state is an n-qubit state, rather than the vac- 
uum state. The new approach requires us to express the 
usual creation and annihilation operators, which act on 
the n-qubit ground state, as a function of the field op- 
erators, which act on the field vacuum state. To derive 
such an expression, we use an explicit model of a single- 



photon source [5] to treat the deterministic creation of a 
qubit from the vacuum. 

The paper is arranged in the following way: in section 
II we show how a qubit algebra can be constructed from 
the vacuum using the creation and annihilation operators 
of a bosonic field with explicit space-time parameteriza- 
tion. In keeping with the Heisenberg Picture throughout, 
we find it necessary in section III to introduce an opera- 
tor corresponding to the deterministic creation of a sin- 
gle particle mode. This allows us to represent arbitrary 
quantum gates acting on a collection of qubits as oper- 
ators on the vacuum state. In section IV we show that 
the algebra is sufficiently general to cope with operations 
on qubits whose wavepackets are not perfectly matched 
to one another. Finally, in section V we illustrate these 
methods by applying them to the simple problem of a 
c-not gate in flat space-time. We conclude by discussing 
the potential applications of this work to problems in 
relativistic quantum information theory. 



II. FIELD QUBITS 

We will consider a field of scalar bosons, which may be 
expanded in terms of its positive and negative frequency 
parts as: 

${x,t) = f dk[G{k,x)a k + G*{k,x)al] 

J (1) 

= A G + Al, 

where G{k, x) is a normalized solution of the Klein- 
Gordon equation and x, k are 4-vectors, with kx = 
g^ v k^x v . We will find it useful to consider wavepack- 
ets, which are localized superpositions of plane wave so- 
lutions: 
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G(k,x) = g(k)e 



ikx 



(2) 



where g(k) is centered on some positive wave number 
kg and is required to be zero for k < 0. Normalization 
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requires: 



/ 



dk | 5 (k)| 2 = l. 



We then interpret the operator Aq as the annihilation 
operator for a particle in the mode G, and A G as the 
corresponding creation operator. From the bosonic com- 
mutation relation [a(k), a'(k')] = S(k — k') we obtain the 
same-time wavepacket commutator: 



i G (fc,x 1; t),i^(/c',x 2 ,i) 
= J dkG(k,xi,t)H*(k,x 2 , 
dk h*(k)g(k)e Xl - X2 . 



t) 



(3) 



In the rest of this section, we will restrict ourselves to 
a set of modes G; that are either perfectly matched in 
their spectral and spatial degrees of freedom, or perfectly 
orthogonal to each other, i.e. 

[A Gi ,A Gj ] = J dk 5i (k) 5 *(k) = %. 

Hence the particular spatial dependence of the modes 
plays no role for the time being, so we will drop the 
subscript G and merely consider the index i, hence 
Ad — > Ai. The spatial dependence will be re-introduced 
in section IV, where we show that a qubit algebra can still 
be constructed in a straightforward manner using infor- 
mation about the particles' spatial profiles; for the mo- 
ment, we describe the construction of a simplified qubit 
algebra from ideal modes with trivial space-time depen- 
dence. To consider qubits, it is convenient to work in the 
Fock basis and write a single particle state as the action 
of A] on the vacuum state: 



II, 



410) 



To build up a qubit algebra, we need to start from a 
'qubit ground state' which is an n particle state in the 
field representation. In the remainder of this section, we 
construct an algebra assuming that we have such an n 
particle state. Section III will then illustrate how parti- 
cle creation can be incorporated in the algebra. Suppose 
that we have n particles, each confined to a pair of or- 
thogonal modes, A\,Bj, where i = 1, 2, .., n labels the ith 
particle. Provided each particle remains confined to its 
two modes, we may interpret the particles as qubits in 
a dual-rail encoding and the state as an n-qubit system. 
For the ith qubit in the two- mode number basis |nA^s)i, 
we define the computational basis states |0), |1): 

|0)i = it|00) = |10) i; 11)4 = ^100)^101)*. 
An arbitrary superposition state is: 

|^) = a|0) i + J 8|l) i = a|10) < +^|01) i ; 



(M 2 + |/3| 2 = l). 
We consider the following operations on the ith qubit: 
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(4) 

These are the quantum Stokes operators [6] . They pre- 
serve the qubit structure of the state space, and within 
it they form a representation of the familiar Pauli al- 
gebra when expectation values are taken. For example, 
even though the operator XX = A^BA^B + B^AA^B + 
A^BB^A + B^ Aw A acting on the ith qubit looks very 
different to / = A* A + W B 1 a quick calculation shows 
that these operators are equivalent under expectation 
values taken in the state | ipi) . It is therefore convenient to 
replace the operators in Q with their equivalence classes, 
which are formally representible by the Pauli matrices: 

1 
1 
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This allows us to perform calculations using just the Pauli 
operators, without having to refer to the fields them- 
selves. At the end of the calculation, we need a recipe to 
convert the final matrix back into field operators in order 
to calculate the expectation value. Motivated by the cor- 
respondence between the field expressions ([4J and their 
matrix counterparts, we propose that a general single 
qubit matrix acting on the ith qubit can be decomposed 
into fields as: 
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P ° 



1 A?A l + SA^Bi + pBi A\ + aB-B t . (5) 



Using the matrix tensor product, we may generalize this 
decomposition to 2™ x 2 n matrices acting on the entire 
rt-qubit system: 

7 5 . . 



P c 



- jA^ A 1 A 2 f A 2 ...Aj A n 
8Ai A^J A 2 ...Aj B n 

... + kBi B\B 2 B 2 --.bJ B n . 



(6) 
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The n-qubit identity in n-qubit space is therefore defined 
by: 



1 



In- 



(7) 



In order to model general n-qubit gates we require a uni- 
versal gate set. We can construct arbitrary single qubit 
operations from the Pauli operators by (for example) 
combining rotations around the Z axis, Z(6) — cosOI + 
isinOZ, with the Hadamard gate, H = l/y/2(Z+X). We 
also need to add a two-qubit controlled sign (c-sign) gate 
to complete the set pQ. In 2-qubit space, the c-sign has 
the matrix representation (in the computational basis): 



U, 



1 



csign 



= - hl 2 + hZ 2 + Zxh - Z X Z<< 



1 

1 

1 





(8) 



In terms of fields, using ([6| and Q, this matrix corre- 
sponds to the operator: 



1 — 2Bi B\Bi B2 



(9) 



Taking this as our c-sign gate completes the universal 
gate set. Arbitrary n-qubit gates can be represented as 
tensor products of the Pauli gates Q and c-sign gate 
([9]) acting on one- and two-qubit subspaces. We can also 
translate between the familiar matrix representations and 
the field operator representations using 



III. SINGLE PARTICLE PRODUCTION 



The qubit ground-state in the preceding section is an 
n-particle state (where n is the number of qubits) and 
so is not the field ground-state. In order to have "pure" 
Heisenberg evolution of the field state we need to identify 
the Heisenberg transformation representing single parti- 
cle creation. Strictly this means finding a unitary trans- 
formation that takes us from the vacuum to the n-qubit 
state. Ideally, we would like to have a unitary Up that 
takes a vacuum mode A as input and gives a single par- 
ticle in the mode Ap as output. This is shown schemat- 
ically in Fig. [l] 
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FIG. 1. Schematic of a unitary single particle source. We will 
find it useful to make use of a large number of ancilla states, 
which begin and end in the vacuum state. 



A. Heisenberg Picture 

We now briefly review how quantum circuit calcula- 
tions can be performed in the Heisenberg Picture [3] . 
Instead of starting with some initial state, \tp) that is 
evolved to some final state, U\i/j) as in the Schrddinger 
Picture, we start from some set of observables rep- 
resenting the measurements made on the system, J, 
and evolve them back to their initial values, U'JU, 
where expectation values can be taken against the initial 
state, {i()\U' JU\i/j) . The measurement operators are con- 
structed from the Pauli operators, i.e. J = f(I, X, Y, Z) 
and, given that the evolution operators are also con- 
structed from Pauli's, i.e. U — g(I, X,Y,Z), calculation 
of U'JU is straightforward from the Pauli commutation 
relations. 

Clearly the two pictures are equivalent with respect to 
expectation values as (^(UUU)^) = {{tp\U^)J{U\^)). 
Here we are interested in "pure" Heisenberg evolution 
for which all evolution, including state preparation, is 
performed in the Heisenberg Picture and so expectation 
values are always taken with respect to the qubit ground- 
state, i.e. (J) = (0\WJU\0). 



Since the input and the output are both pure states, 
such a unitary must exist in principle. We could then 
obtain the output mode as a function of the input 
vacuum modes by calculating its Heisenberg evolution: 

U^AU = A P . (10) 

However, finding an explicit form for the unitary is a 
nontrivial matter because it is a highly non-linear func- 
tion of the field operators. Rather than attempt such a 
task, we will consider a physical model that produces the 
exact same evolution using detection and feed-forwards. 
The model will give us an explicit recipe for writing down 
the output mode Ap in terms of vacuum input modes. 
Specifically, we refer to a suggestion by Migdall, Bran- 
ning and Castelletto [5] for a deterministic single pho- 
ton source that makes use of N spontaneous parametric 
down-converters (SPDCs) and post-selection (see Fig. [2]) 
to produce the desired evolution in the limit of arbitrar- 
ily large N. The basic model is not specific to optics and 
we will present a slightly more abstract version that is 
well suited to our purposes. 
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FIG. 2. A non-unitary near-deterministic single particle 
source. 



A. The Migdall-Branning-Castelletto Source 

Consider the parametric amplification unitary: 
Up = exp { { dkdk'[xP(/fc, k')u\k)v\k') - H.c] 



This unitary acts on the vacuum to create a two particle 
state with probability x 2 ; one particle in the mode u 
and one in the mode v (the u particle may be detected 
to herald the presence of the v particle), x is assumed 
•C 1 and so terms of higher order than \ 2 are neglected. 
The particles have a joint spectral amplitude denoted 
P(k, k') which in general includes entanglement between 
the particle pairs; however for simplicity we will assume 
that the state is separable (achieved in the lab by phase- 
matching) and has the form P(k,k') = G(k)G(k') [7j. 
Under this assumption there are no multimode effects, 
so we can use the shorthand 



dkG(k)e^ kx -^u(k) ^ u , 



and similarly for the mode v. The Heisenberg evolution 
of the operators is then: 



vl = UpuUp = sinh(x)u < + cosh(x)v 
v' = UpvU p = sinh(x)v^ + cosh(x)u . 



(11) 



Generalizing to an array of N such unitaries, labelled 
from 1 to N, we take the i t h unitary to act on the modes 
Ui,Vi with an output amplitude G*(k)G*(k'). The modes 



are assumed to be mutually orthogonal and satisfy the 
commutation relations: 



The output of the array is: 



[v%,v 



Hi 



U Pl ®Up 2 ® ...U P „\0). 



(12) 



Although the probability of any individual SPDC firing 
remains of order % 2 , where x 2 is small, the probability of 
at least one of the SPDCs in the array firing scales with 
N, so that for N sufficiently large, we are almost cer- 
tain to have at least one particle available at any given 
moment. Furthermore, g( 2 > does not scale with N be- 
cause it represents the relative probability of obtaining 
more than one photon for a given mode, independently 
of modes other than the one being considered. Hence the 
probability of obtaining two or more particles in the same 
mode remains negligible. Finally, we perform detection 
on the Ui modes in ( 12 1 and feed forward the topmost 
heralded particle in the array into the mode Ap. All 
other particles that may have been produced are dumped 
and replaced with empty modes. The information about 
the detection result for the mode Hi is contained in the 
particle-number operator: 

n u ,=uJ Uj 

and has eigenvalues 0,1 or 2. Higher values than 2 corre- 
spond to higher orders of x and can be neglected provided 
we choose x to be small, as we do in this paper. In terms 
of the feed-forward and the vacuum modes, we find that 
the output mode Ap has the form: 



N 

£■ 



3-1 
i=0 



1 



cA, 



(13) 



where A is the vacuum input mode referred to at the 
beginning of the section. The expression depends explic- 
itly on the feed-forward n u >, through a bucket-detector 

operator d u i : 

1 



Given the non-negligible eigenvalues of n u 



(14) 

the opera- 
tor ( 14 1 represents a detector that registers a single click 
when one or more particles are present and does noth- 
ing otherwise. This detector is unable to resolve particle 
numbers greater than one; this is in keeping with the 
bucket detector model. 

Any single particle mode must have a particle-number 
expectation value of 1, a negligible probability of con- 
taining more than one photon, and it should satisfy the 
usual commutation relation. Formally, these properties 
are given by the relations: 



(o\ApA P \o) 



1. 



(Q\ApA P ApA P \o) = 0, 
[A P ,Ap] = s(k-k'). 



(15) 
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For Ap to satisfy the commutation relation in (15 1, the 
operator c is taken to be: 

N j-1 
j=l i=0 

We need to verify that Ap also satisfies the other prop- 



erties in (15). From a direct calculation, the details for 



which are shown in the Appendix, we find: 

(4-4|x| 2 + 9|x| 4 )((l-|x| 2 ) Ar ~l) (16) 
5| X | 4 -4 



(h Ap ) = (0\A P A P \ 



9 y ' = 



(o\A P ApApA P \o) 



2|x| 2 (4-5|x 



4\2 



(17) 



(4-4|x| 2 + 9|x| 4 ) 2 ((l-|xl 2 ) 



2W 



In the limit of N large and x small, one finds that 
these quantities approach 1 and respectively (see Figj3|; 
hence the properties ( 15 ) are indeed satisfied. In idealized 



situations, the full expression ( 13 1 is not needed; the iden- 



tities (151 suffice for performing calculations. In more 



practical situations where particle production might be 
non-ideal, or when we may wish to consider relativistic 
observers (see Section VI) , it will be necessary to have 
the explicit expression (13) for Ap in terms of the initial 
vacuum modes. 

Explicitly, if we wish to enact a circuit U G on n qubits 
produced by the above method, where the circuit acts 
on the modes A Gl , A G2 , A Gn , B Gl , B G . 2 , B Gn , then 
instead of calculating the usual general expectation value: 

(^„|F(i Gl , A G2 , ...,B Gl , , (18) 

where \ipin) is the initial n-particle prepared state 

\,l, in ) = A G A G2 ...A G JQ) , 

we now calculate the expectation value: 



(0\F(Ap, Gl) Ap tGa) ...,Bp tGl ,...)\0) 
z (Q\f(u Gl ,u G2 ,...,v Gl ,...)\0}, 



(19) 



where the Ap iGi , Bp Gi are defined in terms of the vac- 
uum modes u Gi , v Gi according to (13). 



IV. PARTICLE MISMATCH 

We would now like to relax the assumption that parti- 
cles in non-orthogonal modes have exactly the same mode 
distribution G(k,x). In the context of wavepackets, this 
will introduce nontrivial commutators of the form ([3]). 
This more general situation might include interactions 
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FIG. 3. (Colour online) Graph of n Ap and # (2) for xi = 0.01 
(red) and \2 = 0.001 (green) as a function of N. The hori- 
zontal scaling has been adjusted relative to xi — 0.01 to bring 
the graphs into coincidence for ease of viewing. The vertical 
scaling is the same for all graphs to allow for a comparison. 
While the graphs of ua p always saturate at tia p — 1 for some 
N, the graphs of </ 2 ' saturate at values that become closer to 
zero (minus infinity on the Log scale) as we make x smaller. 
Thus we can choose x so as to make g' 2 ' as small as we like, 
while still obtaining ua p — 1 for sufficiently large N. 



between particles whose spectral amplitudes are not the 
same, as in quantum optics, or situations in which the 
particles arrive at the gate at slightly different times, or 
their wavepackets only partially overlap in space. In the 
discussion that follows, we will focus on the latter case 
of mismatch caused by mis-timing or spatial discrepan- 
cies, but our analysis applies equally well to other types 
of mismatch. 

One might worry that the algebra of section II no 
longer applies in this scenario, since the extra degrees 
of freedom might allow qubits to become 'invisible' with 
respect to each other, making it unclear how to define 
an n-qubit space. We will make use of a proposal, made 
in the context of optical quantum computation, by Ro- 
hde, Mauerer and Silberhorn (RMS) [5], which allows us 
to separate the amplitudes G(k,x) into parts which are 
separately matched or orthogonal to the interaction of 
interest. We will show that this allows us to retain the 
Pauli algebra in general situations, at the cost of expand- 
ing the state space to include the orthogonal degrees of 
freedom. 

Following RMS, the mismatch between particle modes 
is treated as a rotation into the orthogonal degrees of 
freedom with respect to some 'reference' mode, with the 
amount of rotation being equivalent to the overlap be- 
tween the amplitudes of the modes of interest with the 
reference mode. In this way we can write down an effec- 
tive circuit that can be treated with all the usual meth- 
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ods, yet with sufficient machinery to deal with potential 
qubit mismatch. We will demonstrate this technique for a 
general two-qubit gate. The circuit in Fig|4]is a two-qubit 
system, constructed from two particles with spectral am- 
plitudes G(k,xi) and G(k,X2) and distributed amongst 
the pairs of modes a i ^,b i ^ 1 with i — 1,2 labeling the 

qubits. Hence: Ax,Gi = J dk G(k,xi) d 1 £, etc. The 
qubits interact via a general two-qubit unitary Uc whose 
Hamiltonian is assumed to depend on just a single am- 
plitude H (fc, x c ) representing, for example, an external 
classical pump driving the interaction. Larger circuits 
containing interactions with different pump amplitudes 
can be constructed by connecting basic circuits of the 
type considered here, although each part would have to 
be treated separately. 
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FIG. 4. A general unmatched two-qubit circuit. Here Gi,2 
denotes G(k, £1,2) and H denotes H(k,x c ). 
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FIG. 5. The circuit as seen by orthogonal qubits. 



Our task is to calculate the Heisenberg evolution 
through the circuit given that G(fc, xi), G(k, X2), H{k, x c ) 
are all different in general. We proceed by choosing the 
pump amplitude H(k,x c ) as a reference, and describing 
all other modes by their overlap with H(k,x c ). To this 
end, we will make use of a 'mismatch factor' £, which 
we interpret as being the amount of overlap between the 
interacting wavepackets. Specifically, £ is determined by 
the commutator ^ between the appropriate wavepacket 
creation and annihilation operators at the time of inter- 
action: 



I J dk G(k,x u t)H*(k,: 



(20) 



It turns out that general situations in which neither 
qubit is matched to the reference amplitude can always 
be re-interpeted such that all of the mismatch is car- 
ried by only one of the qubits. For example, if the 
overlap terms are £1 = | J dkG(k,Xi)H*(k,x c )\ 2 and 
£2 = |/ dkG(k,X2)H*(k,x c )\ 2 , we can multiply them 
to form a third quantity £3 = £1^2 and then interpret 
£3 as the mismatch of just one qubit whose partner is 
perfectly matched to the reference amplitude. Hence we 
lose no generality by assuming that one of the qubits is 



perfectly matched and the other is not. We will there- 
fore let G(k,X2) — H(k,x c ) and let the entire mismatch 
be carried by the overlap of G(k, Xi) with H(k,x c ). We 
first write the modes as superpositions of matched and 
orthogonal components: 



i ljGl ee y/£Ai tH + Vl-£c, 
Bi, Gl = \flK11 + Vl^d, 



(21) 



where we have introduced the auxiliary vacuum modes 
c, d. From a circuit perspective, we may view the G\ 
modes as the result of a rotation that causes mixing be- 
tween the matched and orthogonal parts. Formally, this 
is achieved by the evolutions: 



c \fic - y/l-£ax, 
h -> VI h + d, 
^d- y/l^ih, 



(22) 



where it is understood that the modes cij, bi are matched 
to H{k,x), i.e. we identify dj = A^H{k,x), etc. The 
choice of phase convention in the above evolution doesn't 
matter because the modes representing the unmatched 
part are initially empty. We note that even if the spectral 
profiles of the interacting modes are identical, h(k) = 
g(k), spatio-temporal differences between the modes such 
as that produced by a delay on one qubit may still give 
rise to a non-trivial overlap, since G(k, x±) ^ G(k, x%) in 
general. We recover the algebra of section II as a special 
case: when £ = 1, then qubit 1 is perfectly matched 
and the orthogonal subspace is trivial; conversely if £ = 
then qubit 1 is completely orthogonal and the qubit 
subspaces decouple. 

The gate Uc only acts on the matched modes 
Si, 61, 02, 62- The orthogonal modes c, d encounter a dif- 
ferent circuit, U c , that can be obtained from Uc by phys- 
ical considerations, as we now show. To this end, we 
make use of the fact that any two-qubit unitary can be de- 
composed into single qubit gates and c-signs. Specifically, 
the general two-qubit unitary Uc can be decomposed 
into eight single-qubit operations and three c-signs [5], 
as shown in Fi g|4| The single-qubit gates may be con- 
structed from (41) as normal, using field modes matched 
to G{k,x\) and H(k,x c ) for qubits 1 and 2 respectively. 
As the single-qubit operations are independent of one 
another, it makes no difference whether the qubits are 
matched or not; commutators between unmatched modes 
simply do not arise. The two-qubit c-sign gates are more 
complicated; we must make reference to a particular 
physical model. A natural choice, applicable to optical, 
microwave and ionic aubits flOl II lj . is the strong nonlin- 
ear cross-Kerr effect between two modes. To the external 
pump mediating this effect we prescribe the amplitude 
H(k). The unitary operator for the c-sign acting on the 
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modes bi,b 2 is: 



UKerr = exp[-iirh bl h b2 ] 

= 1 + [e~™ - l]n bl n b2 
= 1 - 2n bl n b2 

= 1 — 2_B^ 1 Bh,iB\ { 2 Bh.2 j 



(23) 



where we have used the property: 

(n bl h b2 ) p = n bl n b2 , 

which follows from that fact that the modes b\ , b 2 can 
each only contain or 1 particle. Note that the unitary 
( 23 ) is consistent with an abstract c-sign ^ whose modes 
are matched to H(k,x c ). From Q we see that the ab- 
sence of one or both particles implies U CS i gn — > 1. This 
implies that the circuit U' c seen by the orthogonal modes 
c, d is just the decomposition of Uc without the c-signs, 
see Fig{5] We can now interpret the evolution (22) as the 
action of a gate U r on the joint system Uc 
will be useful as we now show by example. 



U^, which 



V. AN EXAMPLE: THE C-NOT 

To illustrate the use of the techniques described above, 
we will calculate the output of a simple two-qubit cir- 
cuit in which the qubit pulses have identical spectral 
properties but do not fully overlap at the gate. This 
could be achieved by placing a slight delay on one of 
the pulses, leading to a spatial discrepancy between the 
propagated modes at the time of interaction. We con- 
sider a circuit consisting of a single controlled NOT 
(c-not) gate, see Fig{6j The first qubit is prepared 
in an arbitrary superposition by applying the operator 
U s = yL — o?l\ — cdY\. The output of this circuit in the 
Schrodinger picture for the ideal matched case is the Bell 
state: \/l — a 2 |00) We will now investigate how 

the spatial mismatch affects the functioning of the gate. 
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FIG. 6. A c-not between unmatched qubits, with an arbitrary 
rotation on qubit 1. The c-not gate is decomposed as a c-sign 
and two Hadamard gates. 



For this example we will assume that the qubits are 
spectrally matched to each other and the gate; g(k) = 
h(k). The mismatch is introduced by defining G(k, X\) = 
.g(fc)e (4fel ~ wt) for qubit x and G{k,x 2 ) = g^e^-^ 
for qubit 2 and the gate interaction, where the modes 
are now the same apart from being centered at different 
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FIG. 7. The circuit as seen by orthogonal qubits. 



points in space, X\ ^ x 2 , 
that the mismatch £ is: 



at time t. From (20), we find 



dkg{k)g*(k)e l ^ Xl 



-x 2 )|2 



(24) 



The degrees of freedom are denoted by the orthogonal 
modes: 



{fii, 6i,c,d, 02,62} • 



(25) 



The complete circuit containing U cno t is shown in Fig|6j 
whereas the circuit seen by the orthogonal modes is 
shown in FigjTJ The c-not gate in the latter circuit, which 
we denote U' cnot , is found to be the two-qubit identity 
using the methods of the preceding section. The modes 
ai,&2 are single particle modes and have the form (13) 
because the initial state is |00), which corresponds to the 
creation of a particle in each of the modes a\ and 0,2- The 
detectors are taken to be bucket detectors. Hence a mea- 
surement of the expectation value of the Pauli operator 
X on qubit 1 is represented by calculating the expecta- 



tion value of the operator + 



x[ cd ^ where 



etc. 



All unitary operators acting on the qubit 1 subspace can 
be similarly decomposed into a sum of matched and or- 
thogonal parts, corresponding to the ai,b\ and c, d sub- 
spaces respectively. Given the algebra (Ek together with 
knowledge of the orthogonal circuit Fig(7J we can con- 
struct a matrix representation of general two-qubit op- 
erations in the space spanned by {d\, b\, c, 0?} ® {02, 62} 
using 8x8 matrices in block diagonal form in which the 
two blocks represent the matched and orthogonal sub- 
spaces: 



U c 



U c 
U' c 



To avoid confusion, we will denote operators acting on 
the entire space with a bar on top, J\, instead of a hat, 
Ji , reserving the latter for operators acting on either sub- 
space (since qubit 2 is perfectly matched, the unmatched 
subspace is trivial for that qubit: K 2 — K 2 )- We now 
write down the matrices for the gates in the circuit ([6]): 



U rnn t 



Ucnot 

/1/2 



(where U s 



U S I 2 
UJ 2 



\f\ — a 2 —a 



a 



vT 



or 



), 



U r 



sflhU 



Intuitively, the matrix forms of U cno t and U s can be un- 
derstood in terms of their action on the matched and 
orthogonal subspaces; the c-not acts differently on each 
subspace due to the varying behavior of the non-linear C- 
SIGN gate, but the rotation is just a linear single-qubit 
gate and works regardless of whether the qubit has been 
rotated into orthogonal modes or not. The unitary U r 
is simply the rotation of the Si, 61 modes into the c,d 
modes according to ( 22 1 . We can now calculate the ex- 



pectation value of some arbitrary operator J on qubit 1 
and K on qubit two. This operator may be represented 
by the matrix: 



•hK 2 



JiK 2 
J X K 2 



We wish to calculate the quantity 

(OlUU^UsyslO) 



(26) 



where U sys is the combined evolution of the circuit, 



(27) 



We find: 



U S y S 



V£,U cno tUib s I 2 — £,U C notU\bsI 2 



For a given J and K 1 we can then obtain the matrix: 



(28) 



Taking the expectation value of this operator in the vac- 
uum is trivial since the only modes containing particles 
are the 64 and a 2 modes, so all matrix elements contain- 
ing other terms annihilate the vacuum and vanish. The 
only remaining term is the coefficient of 



a}aia 2 ci 2 |0) = 1, (using the identities (15)) 



whic h co rresponds to the top left hand entry in the ma- 
trix ( 28 1 . In summary, once we know U sv s , w e find the 



expectation value of J\K 2 by calculating (28) and then 



reading off the top left corner entry. For the case where 
J 1 = I\, K 2 = Z 2 we obtain the result: 

(hZ 2 ) = l-2a 2 £. 



Depending on whether £ is closer to zero or one, this term 
approaches either 

(o 1 o 2 |/ 1 z 2 |o 1 o 2 ) = 1, 



or 



(0 1 2 |Z 1 Z 2 |0 1 2 ) = I -2a 2 . 

These two limits correspond to the c-not either not work- 
ing or working, depending on whether the qubits meet in 
the gate; that is, on how well their individual amplitudes 
overlap during the interaction. In this example, accord- 
ing to ( 24 ) , £ depends on the shift in position xi — x 2 of 



the pulses. 



The key novelty of our formalism is the expression ( 28 1 , 



which incorporates the dynamics and spectral features of 
the qubits into the operator algebra, and the fact that 
these expressions may be written out explicitly in terms 
of field operators acting on the global vacuum according 
to the procedure outlined in section III (i.e. the expres- 
sion has the form (19)). 



VI. CONCLUSIONS AND FUTURE WORK 

We have developed tools for doing quantum circuit cal- 
culations using the particles of a bosonic quantum field 
for qubits. The algebraic structure of the circuit is de- 
fined with respect to a global vacuum state, and expecta- 
tion values are taken in this state in accordance with the 
Heisenberg Picture. The usual techniques for calculating 
evolution on quantum circuits, (matrix representation or 
Pauli operator algebra) apply here, but now include the 
method of RMS [8] for dealing with particle mismatch. 

A strong motivation for our work comes from the new 
area of relativistic quantum information theory (RQI) 
- concerned with adapting the tools of quantum infor- 
mation to relativistic settings (eg. [T2Tl20j ) . Typically, 
observers in different (especially non-inertial) reference 
frames can observe inequivalent ground states for fields 
shared between them [HI [2TJ [22]. Calculations for re- 
alistic scenarios can rapidly become intractable in the 
Schrodinger Picture due to the highly non-trivial trans- 
formations of the ground-state seen by different observers 
that may be required. In contrast, working in the Heisen- 
berg Picture means the ground-state can be chosen in a 
suitable inertial frame and all the transformations be- 
tween frames can be achieved via more straightforward 
Bogolyubov Transformations of the operators 23 . 

More generally it is important that formulations of 
quantum information problems are consistent with cur- 
rent methods for doing quantum field theory on curved 
spacetimes, particularly the algebraic approach to quan- 
tum field theory [24l I25j. This approach is based on an 
algebra of observables constructed from products of field 
operators, which can take on the interpretation of particle 
creation and annihilation operators acting on a vacuum 



9 



state. Again, the formalism is expressed in the Heisen- 
berg Picture and expectation values are taken in the ap- 
propriate vacuum state. The formalism we describe here 
can be considered as a first step towards developing a 
generalized qubit algebra that is able to incorporate gen- 
eral transformations in the relativistic setting. 
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Appendix A: Calculation of (fiAp) and 

In this appendix we give the details for calcul atin g the 
expressions for (fiA P ) and g*- 2 ' in Eqs. ( |l6[ ) and (17) from 
Section |III A| We consider the general SPDC evolution 
given in Eq. ( 11 ), using the shorthand notation 



r.t 



lj = pvj + quj 



fit 



Uj = puj + qv'j 



(Al) 



instead. When calculating (0|A P Ap|0) and 

(0|Apv4p74pAp|0) we consider a more general de- 
tector model of the form 



d u >. = (pi - vh u >)h v 



(A2) 



which is equivalent to Eq. (14) when /i = 3/2 and v 
1/2. 



1. Calculation of (0\A P A P \0) 

We write the full expression for (0|A pJ 4p|0) given Ap 
in Eq. (pi: 



N fc-1 j-l 

(EI!' 1 - d u ,)v'^d u , k d u , fij II (1 - d< )> , (A3) 

j,k=l 1=0 i=0 

where it is understood that all expectation values are 
taken with respect to the vacuum. The summation can 
be broken up into three separate terms: 



N N N N 

E = E + E + E 

j,fc=l 3<k=l ]=k=l j>k=l 



(A4) 



Since {v'^d u ' k ) = 0, only the J2j=k term in Eq. (A3) 



survives. Given that 



(vjd u >. d u >. iij) = fi (p, q, fi, v) 
((1 - d u >) 2 ) = f 2 {p,q,H,v), 



(A5) 



where fi(p, q, fi, v) and f2(p, q, v) are both polynomial 
functions of p, q, \x and is, we find: 



N j-1 

(AUk)A P (k)) = h E II h = h 

j=l 1=0 



(A6) 



When we substitute \x = 3/2, v = 1/2, p = 1 and q = \x\, 
where |x| <C 1 for a SPDC, we obtain Eq. (16) to order 



\X\ 



2. Calculation of {0\A P A P A P A P \0) 

We write the full expression for (ApApApAp) given 

Ap in Eq. (13) using the fact that [v^ , d u >] = [v^d-u'] — 
0: 1— r ' 



JV 



fc-i 



E d k dw k «a du' a d u r v'j d K v'cYii 1 - d u ' t ) 



•=i 



1=0 



j-1 



c-1 



(A7) 



b=0 



d=0 



where again it is understood that all expectation values 
are taken with respect to the vacuum. The summation in 
this case can be broken up into 75 separate terms, most 
of which do not contribute: 

• The 24 permutations of summations of the form 
J2 a >c>j>k d° n °t contribute since {v k d u ' k ) = 0. 

• The 12 x 3 = 36 permutations of summations of the 
form £ 0=c>J>fe , £ >c=j>fc and E a >c>j=fc do not 
contribute since (v^d u ' v^d u i ) — and (v' k ]d u > (1 — 
^) 2 )=0. 

• The 4x2 = 8 permutations of summations of the 
form Y,a=c=j>k and E a >c=i=fc do not contribute 
since {v£ d u ' k v£ d u > k d u > k v' k ) = 0. 

• Four of the 6 permutations of summations of the 
form X) a =c>j=fc are n on-zero and given by 

a-l fe-i 

E (^d u ^d u ,j u ,j'j u , k v' k - d u ,f j](i - d u[ f) 

a=c>j=k b=0 1=0 

fc-1 a-1 

= E^II/3 II /a, (A8) 

a>k 1=0 b=k+l 



where 



((1 - d U ' k ) A ) = hip,q,n,v) 
(v'f f .d Uk d u ' k v' k (l - d ul J) = f 4 (p,q,n,iy) 

are also polynomial functions of p, q, \x and v. 



(A9) 
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• The summation X)a=c=j=fc ^ s non-zero and given 

by 

fc-i 

a—c—j—k 1—0 
k-l 

=E/^n^ ( ai °) 



k 1=0 



is a polynomial functions of p, q, fi and v. 
The expression for (ApApApAp) therefore reduces to: 



k-l a-x 



N k-l 



(A12) 



a>k 1=0 b=k+l k=l 1=0 

= 4/l/4 



h-h + f?{h 1) - /f(/3 1 ^ . /; /7 3 * - 1 



(/a - l)(/ 3 - l)(/a - h) 



/3-1 



where 



When we substitute = 3/2, f = 1/2, p = 1 and q = \\\, 
where \x\ < 1 for a SPDC, we obtain Eq. (17) to order 
(Vkdu'Jkdu'Ju> k v k d u > k v' k ) = f 5 {p,q,n,v) (All) |x| 4 . 
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